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ABSTRACT 

High-resolution ultraviolet observations of the black hole X-ray binary Cygnus 
X-1 were obtained using the Space Telescope Imaging Spectrograph on the Hub- 
ble Space Telescope. Observations were taken at two epochs roughly one year 
apart; orbital phase ranges around 0orfe = and 0.5 were covered at each epoch. 
We detect P Cygni line features from high (N V, C IV, Si IV) and absorption lines 
from low (Si II, C II) ionization state material. We analyze the characteristics 
of a selection of P Cygni profiles and note, in particular, a strong dependence 
on orbital phase for the high ionization material: the profiles show strong, broad 
absorption components when the X-ray source is behind the companion star and 
noticeably weaker absorption when the X-ray source is between us and the com- 
panion star. 

We fit the P Cygni profiles using the Sobolev with Exact Integration method 
applied to a spherically symmetric stellar wind subject to X-ray photoionization 
from the black hole. Of the wind-formed lines, the Si IV doublet provides the 
most reliable estimates of the parameters of the wind and X-ray illumination. The 



^ Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorpo- 
rated, under NASA contract NAS5-26555. These observations are associated with programs GO-9646 and 
GO-9840. 
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velocity v increases with radius r (normalized to the stellar radius) according to 
1^ = i^ooi^ ~ fi^l^Yi with j3 ~ 0.75 and v^o ~ 1420 km s~^. The microturbulent 
velocity was ~ 160 km s~^. Our fit implies a ratio of X-ray luminosity (in units 
of 10^^ erg s"^) to wind mass-loss rate (in units of 10~^Mq yr~^) of Lx.ss/^-e ~ 
0.33, measured at M_6 = 4.8. The lines from the lower ionization species and 
the He II A 1640 absorption are consistent with formation in the photosphere of 
the normal companion. 

Our models determine parameters that may be used to estimate the accretion 
rate onto the black hole and independently predict the X-ray luminosity. Our pre- 
dicted Lj; matches that determined by contemporaneous RXTE ASM remarkably 
well, but is a factor of 3 lower than the rate according to Bondi-Hoyle-Littleton 
spherical wind accretion. We suggest that some of the energy of accretion may 
go into powering a jet. 

Subject headings: black holes binaries: individual (Cyg X-1) — ultraviolet: stars — 
X-rays: winds 

1. Introduction 

Cygnus X-1, the only Galactic X-ray binary with a high mass companion where existing 
observations require a black hole for the compact object, was first discovered in 1962 (Cowley 
1992, and references therein). In addition to the well-established high mass function found 
with optical observations, the X-ray data of Cyg X-1 display transitions from a high flux 
state in the 2-10 keV band (where a strong soft component dominates) to a low flux state 
(where the soft component largely disappears) that has been interpreted as characteristic of 
black hole systems. It also displays highly broadened Fc Ko; emission (Miller et al. 2002) 
that is consistent with models for X-ray reflection in Galactic black holes and AGNs. The 
broad-line shape of Fe Ka may be caused by Doppler shifts and the gravitational field of the 
black hole. 

The Cyg X-1 system consists of a supergiant star and a compact object. The mass of the 
compact object is in the range 7-20 Mq (Shaposhnikov & Titarchuk 2007; Ziolkowski 2005) 
the mass of the visible star is in the range 18-4OM0 (Ziolkowski 2005; Tarasov et al. 2003; 
Brocksopp et al. 1999). The binary orbital period is 5.6 days (Bolton 1972). Miller et al. (2005) 
using Chandra/HETG observations find that the X-ray spectrum of Cygnus X-1 at phase 
0.76 is dominated by absorption lines, in strong contrast to spectra of other HMXBs such as 
Vela X-1 and Cen X-3. Schulz et al. (2002) report marginal evidence for ionized Fe transi- 
tions with P-Cygni type profiles at orbital phase 0.93 whereas Marshall et al. (2001) find no 
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evidence for such line profiles at phase 0.84. Miller et al. (2005) suggest that, while the spec- 
tra of Cen X-3 can be modelled by a spherically-symmetric wind (Wojdowski et al. 2003), the 
X-ray absorption spectrum of Cyg X-1 requires dense material preferentially along the line 
of sight; considered together, the Chandra spectra provide evidence in X-rays for a focused 
wind in Cygnus X-1. 

Our Space Telescope Imaging Spectrometer (STIS) observations of Cygnus X-1 were 
obtained when Cyg X-1 was in its soft/high X-ray state and show line profiles that change 
significantly between orbital phases 0.0 and 0.5. We interpret these changes in terms of 
models that include the effects of X-ray photoionization on the stellar wind of the normal 
companion or of a focused wind. We test our model predictions with contemporaneous X-ray 
observations. The observations and analysis are described in Section 2, our models of the 
line profiles are presented in Section 3, and our interpretation and conclusions are discussed 
in Section 4. 



2. Observations and Analysis 

Cyg X-1 was observed with the STIS on the Hubble Space Telescope (HST) when the X- 
ray source was behind the normal star and half an orbit later on two separate epochs roughly 
one year apart. Figure [1] shows the times of our observations in comparison with the 1 day- 
averaged light curves obtained from the all-sky monitor on the Rossi X-Ray Timing Explorer 
(RXTE; Levine et al. 1996). The RXTE light curves show that the HST observations were 
taken when the binary was at relatively high X-ray flux which is associated with the X-ray 
soft/high state. 

The STIS instrument design and in-orbit performance have been described by Woodgate 
et al. (1998) and Kimble et al. (1998). The E140M grating provided a resolving power of 6 
km s'^ in the wavelength range 1150-1740A. The data were processed through the standard 
HST/STIS pipeline and further reduced using the STSDAS routines available through IRAF. 

Table 1 lists the HST observation identifier, along with the start dates, exposure length, 
and orbital phase at the start of each observation. Phases were calculated using an ephemeris 
that places phase zero at JD2441874.707±0.009, with an orbital period of 5.599829±0. 000016 
(Brocksopp et al., 1999). Phase zero corresponds to supergiant inferior conjunction. 

Figure [2] shows the eight datasets in increasing orbital phase order regardless of epoch. 
These data are not smoothed and not corrected for reddening and no instrumental quality 
control has been applied (i.e., removal of hot pixels, etc). It is clear that the changes in high 
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ionization material line profiles are orbital phase dependent, with stronger absorption near 
orbital phase 0.0. The close-ups of the N V, Si IV, and C IV regions in Figure 3 demonstrate 
that this orbital phase dependence persists for observations taken a year apart. 

Figure m shows a spectrum taken at orbital phase 0.96, with the stronger spectral lines 
labeled. Table 2 lists the features identified using the program Spec View U and tables sup- 
plied by NIStB We note that some of the lines are saturated and several are blended with 
other features. Figures IMTUl show profiles of both high (N V, C IV, Si IV) and low (Si II, 
C II) ionization state material in each of two orbital phases at two epochs. The low ioniza- 
tion material and He II are likely to be from the stellar surface rather than the wind. 



3. Line Profile Models 

We have implemented a model that uses the Sobolev method with Exact Integration 
(SEI; Lamers et al. 1987) to predict the P Cygni line profiles from a wind ionized by an 
embedded X-ray source. The details of the wind ionization are identical to those given 
by Boroson et al. (1999; equations 11-16). The SEI method extends the escape probabil- 
ity method of McCray et al. (1984) by including integration along the line of sight. This 
allows scattering from a range of points around the "resonant point", as will result when 
there is small-scale turbulence (microturbulence). The SEI method also allows a more exact 
treatment of interaction between P Cygni line doublet components. 

We take note of a complementary analysis of the same STIS observations by Gies et al. (2007). 
Gies et al. also use the SEI method to analyze the wind, but assume that the ions that form 
the P Cygni lines are only present in the region of the wind in which the primary blocks 
X-ray ionization from the black hole. 

Following McCray et al. (1984), we compute the local ionization fraction by combining 
ambient photoionization (from shocks in the wind, etc.) that would be present in an isolated 
O star with the ionization rate computed using the XSTAR code, assuming a local ionization 
parameter ^ = Lx/nr^ (with X-ray luminosity L^, number density n, and distance from the 
black hole r^). We use the X-ray spectrum shape modeled by Wilms et al. (2006) for 
observations taken near our HST observations. Where the O star shadows the wind from 
X-ray photoionization by the black hole, we assume only the ambient photoionization rate. 



^ http : / / www. stsci . edu / resources /softwareJiardware / specview 
'http://www.physics.nist.gov/PhysReiData/Handbook/index.html 
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which gives optical depths in the hne parameterized by ai and a2, or Twind, often used for 
winds of isolated OB stars (Lamers et al. 1987). The XSTAR ionization code is described in 
Bautista & Kallman (2001). 

We fit the model, which includes adjustable free parameters, to the spectra at phases 
near both superior and inferior conjunction of the black hole. Thus, the fits are sensitive 
to both the bhic-shiftcd absorption trough and the red-shifted emission peak. Although the 
lines are saturated, the optical depths of the lines in the absence of X-ray photoionization 
still affect the resultant fits. 

We attribute the sharp features (velocity widths of order 30 km s~^) to interstellar 
absorption lines, the narrow absorption lines (widths of 94 km s~^) originate on the stellar 
photosphere, and the broadest features come from the wind (terminal velocity from our 
model of 1420 km s~^). In addition, when we model the NV lines, we incorporate narrow 
interstellar Mgll lines at 1239.925, 1240. 397A. The SEI calculation of the line profiles near 
the rest velocity can require integration to start extremely close to the stellar surface in order 
to agree with the results of comoving frame methods (Lamers et al. 1987). Hence we expect 
the P Cygni line fits to be least rehable near the rest wavelengths of the lines. 

We assumed fixed values for parameters describing the O star and orbit as given in Table 
3. Our results are not sensitive to the abundance values used, except for the determination 
of M. The parameters describing the wind and X-ray illumination are allowed to vary until 
a minimum is found through the downhill simplex ( "amoeba" ) method of Nelder & Mead 
(1965). The errors we quote are the rms variation in the final simplex, in which the xl 
varies by less than 1%. There may be other local minima and some of the parameters 
may interact. For example, different values of Twind, oci, a2 may produce similar functions of 
optical depth versus radius in the wind. 

The Si IV doublet is generally the most reliable indicator of wind behavior in stars 
of this spectral type (Lamers et al. 1999), as the NV and CIV lines are saturated. We 
show that the lines can be fit using a wind acceleration law that is standard for OB stars: 
V = Voo{^ — with r the radius in the wind normalized to the stellar radius, and v^o and 

/3 free parameters of our model. (We note that our code used v = 0.01^00 + 0.99*^00(1 — 1/r)'^ 
in order to avoid singularities). The best- fit value of /3 is 0.75, which is standard for OB 
star winds. The best-fit value of Voo ~ 1420 km s~^ is lower than the 2300 km s~^ found by 
Davis & Hartmann (1983). It is possible that the higher terminal velocity found by Davis & 
Hartmann was due to line blending on the blue side of the feature that was hard to discern 
in the weak lUE spectrum. We also note that Davis & Hartmann used lUE data of mostly 
C IV to determine their terminal velocity whereas we use STIS data of Si IV. Also we allow 
for microturbulence. For the CIV lines we fix (3 to that found for the Si IV fits, and then for 
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N V, which has a lower signal to noise ratio, we fix also foo- The best-fit parameters for the 
fits shown in Figures [TllfT3l are given in Table 4. 



4. Accretion Rate and X-ray Luminosity 

4.1. Mass Accretion Rate 

The fits to the changing P Cygni lines by a model of X-ray ionization within a spherically 
symmetric wind determine parameters that may be used to estimate the accretion rate onto 
the black hole. For accretion purely through capture of stellar wind material, Bondi & Hoyle 
(1944) show that the rate at which mass is captured from a stellar wind by a compact object 
(^capture) IS giveu by 

-"^-'capture V / 

where G is the gravitational constant, Mbh is the mass of the compact object, p is the 
density of the undisturbed gas fiow near the compact object, and V^rei = ("^wind + '^orbit)^''^ 
the velocity of the wind relative to the compact object, which has velocity t'orbit = 27r_R/P 
for orbital radius R and period P in a wind we take to have fwind = VooiX ~ R*/RY ■ From 
mass conservation we have 

- ^ 

We can relate Mcapture to the X-ray luminosity if we assume accretion releases energy 
with some efficiency e ~ 0.1, so that = eMcaptmeC^ ■ For disk accretion, we expect 
e = 0.057 for a nonrotating black hole and e = 0.42 for a maximally rotating black hole 
(Shapiro & Teukolsky 1983, p. 429). 

Thus our best-fit values for M, /?, and Voo provide an independent prediction of L^. 
This prediction is based on the idealization that the accretion proceeds entirely through 
the gravitational capture of stellar wind material. The prediction also depends on other 
uncertain parameters (e, Mbh). Our estimate, ^ 8 x 10^^ erg s~^ depends sensitively on 
the velocity of the stellar wind near the black hole. As Voo varies from 1000 to 2000 km s~^, 
Lx goes from SxlO^'' to 3 x 10^^ erg s~^. 

The model dependence on the properties of the X-ray ionized region of the wind also 
leads to a determination of L^/M and M, which together give = 1.6 x 10^'' erg s~^ and 
M = 4.8 X IO^^Mq yr-i in the fit to the Si IV lines. The CIV and NV lines, which we 
expect to be less reliable, give L^^ = 1.8 x 10^^ and L^. = 8.0 x 10^"^ erg s~^, respectively. 
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We can compare these models for the X-ray luminosity with the observed contempora- 
neous RXTE ASM count rate, ~ 80 counts s~^ (Figured]). Schulz et al. (2002) observed 
Cygnus X-1 with Chandra during a period in which the RXTE ASM showed flares. An upper 
estimate for the ASM count rate during their observation, 50 cts s~^, together with their 
measure of the 0.5-10 keV X-ray luminosity as 1.6x10^'' erg s~^ (for a distance of 2.5 kpc), 
would imply an X-ray luminosity of > 2.6 x 10'^'' erg s~^ during the STIS observations. 



4.2. X-ray Luminosity 

One further uncertainty in our model of the wind is that we assume a constant X-ray 
luminosity whereas the X-rays from Cygnus X-1 often flare and display a complex power 
spectrum. The saturated P Cygni lines may respond nonlinearly to flares, although light 
travel times may diminish the response. In Figures [H] and [151 we show how the P Cygni 
lines may change in our model in response to a change in X-ray luminosity. If the change 
were linear, the graph of the model with best-fit value of would be identical to the average 
of the other two graphs. Whereas we see a clear asymmetry. 

In Figures [T6|, [TTl and [181 we show contours of constant values of Logio(asiiv)5 where 
asiiv is the fraction of Si that is in the form Si IV. The lines extending radially from the 
center of the O star show the lines of sight at = 0.55 (extending to the bottom) and 
= 0.96 (extending to the top), adjusted for orbital inclination i by the use of an effective 
phase (j)' such that cos27r0' = cos27r0sinz. 

Figures [TH] through [TS] show that the black hole is very effective at removing the Si IV 
ion from the wind. However, it is difficult to see intuitively the optical depth in the wind, 
given the ion fraction. Therefore, we show in Figures [T9] through [2T] the radial optical depths 
of Si IV in the wind, given the same three X-ray luminosities, Lr^ = (1/3, 1, 5/3) x 1.6 x 10^'^ 
erg s~^, with other parameters fixed to the best-fit values. From these plots it should be 
apparent that the global ionization in the wind is sensitive to the X-ray luminosity. This is 
particularly true at the low end of the X-ray luminosity range. In that case, a large region 
outside of the X-ray shadow remains at an optical depth that causes noticable signatures in 
the line profiles, in spite of X-rays from the black hole. 

In these figures, we use our parameterization of the background optical depth in terms 
of ai, a2, and r„ind to determine the optical depth in the shadow region. The resulting 
contours of constant optical depth are ellipsoidal, compacted towards the compact object. 
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5. Discussion and Conclusions 

The Space Telescope Imaging Spectrograph on Hubble provides the highest resolution 
ultraviolet spectra taken of Cyg X-1 to date. Observations were taken at two epochs roughly 
a year apart: at each epoch orbital phases when the compact object is behind the stellar 
companion and when the compact object is in front of the companion star were covered. 
We find P Cygni profiles from high ionization (NV, CIV, Si IV) gas. For both epochs the 
P Cygni profiles show significantly less absorption at phases when the compact object is 
in the line of sight. RXTE observations indicate that the X-ray fliix of the system was at 
a similar level at each epoch. The observed changes can be attributed entirely to orbital 
effects. We interpret this to mean that X-rays from the compact object photoionize the wind 
from the massive companion resulting in reduced absorption by the wind material. 

P Cygni profiles of selected species are consistent with the Hatchett-McCray effect, in 
which X-rays from the compact object photoionize the stellar wind from the companion star, 
thereby reducing absorption. This effect also appears in UV observations of LMC X-4 and 
SMC X-1 (Boroson et al. 1999; Vrtilek et al. 1997; Treves et al. 1980). SEI models can 
fit the observed P Cygni profiles and provide measurements of the stellar wind parameters. 
The Si IV fits are the most reliable and we use them to determine L^^/M = 4.8 ± 0.3 x 10^^ 
ergs s"-*^ Mq-*^ yr, where Lx is the X-ray luminosity and M is the mass-loss rate of the star. 
The results from the C IV and N V hues are less rehable because they are saturated and the 
C IV fit does not match the data well. For these fits we fixed the terminal velocity i^oo and 
the microturbulent velocity to those given by our fits to Si IV. 

The best fit values for the optical depth in the ambient wind are high (> 10). Once 
saturated, the OB star wind lines hardly change with large optical depth, but when much 
of the wind is ionized by the black hole, the ion fraction in the remaining regions can have 
a significant effect on the line profile. 

The reduced absorption when the compact object is in the line of sight is inconsistent 
with focusing of the wind toward the compact object, as has been suggested by several 
authors (e.g. Sowers et al. 1998; Tarasov et al. 2003; Miller et al. 2005), as then we would 
expect more absorbers in the line-of-sight and hence increased P Cygni absorption. Also, 
lUE observations taken at 8 orbital phases show a continuous variation in the P Cygni 
profiles with maximum absorption at phase 0.0 and minimum at 0.5 (Treves et al. 1980; van 
Loon et al. 2001). Further high spectral resolution ultraviolet observations of Cyg X-1 will 
be necessary to study the behavior of the P Cygni fines during different X-ray states. 

In an analysis of 2 years of RXTE/ASM data Wen et al. (1999) found tfie 5.6 day orbital 
period of Cyg X-1 during tfie X-ray low/fiard state, but no evidence of tfie orbital period 
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during the high/soft state. Wen et al. suggest that absorption of X-rays by a stellar wind 
from the companion star can reproduce the observed X-ray orbital modulations in the hard 
state: The lack of modulation in the soft state could be due either to a reduction of the 
wind density during the soft state or to partial covering of a central hard X-ray emitting 
region by an accretion stream. Gies et al. (2003) used the results of a four year spectroscopic 
monitoring program of the Ha emission strength of HDE226868, the normal companion 
to Cyg X-1, to argue that the low/hard X-ray state occurs when there is a strong, fast 
wind and accretion rate is low, while in the high/soft state a weaker, highly ionized wind 
attains only a moderate velocity and the accretion rate increases. The interpretations of 
both Wen et al. (1999) and Gies et al. (2003) are inconsistent with the fact that the total 
X-ray luminosity from 1.3-200 keV remains constant during both the X-ray soft and hard 
states (Wen et al. 1999): the designation of X-ray high or X-ray low during these states is 
only applicable for the 1.5-12 keV ASM band. Since the 1.3-200 keV X-ray luminosity is 
unchanged from the hard to soft state, fluctuations in the narrow ASM band cannot be due 
to reduction in accretion, or obscuration of the X-ray source; rather it is a physical change 
that causes the dominant emission mechnism to switch from thermal to power-law. 

We note that the orbital modulation observed by Wen et al. in the hard state (±1.6 
ASM cts/sec around the average) is less than the errors on the counts during the soft state 
(±3 ASM cts/sec); we suggest that the quahty of the RXTE ASM is not sufficient to detect 
this low level orbital modulation during the soft state. Ultraviolet observations clearly show 
orbital modulation during both X-ray hard and soft states (Gies et al. 2007). Our wind 
models explain both the X-ray and ultraviolet flux during the soft/high state. We need high 
spectral resolution ultraviolet observations during the hard/low state to determine if there 
is a change in wind density between states. 

Our determination of the mass-accretion rate can be considered as a positive check on 
the Hatchett-McCray models, as the Si IV line gave = 1.6 x 10^^ erg s~^, and the model 
is subject to systematic uncertainties. However, this value of is a factor of 3 lower than 
the best estimate of the accretion rate according to Bondi-Hoyle-Littleton spherical wind 
accretion. We suggest that some of the energy of accretion may go into powering the jet. 

Our test of the dependence of our results on X-ray luminosity conflrms the utility of our 
models and demonstrates that the wind outside the shadow zone is still sensitive to X-ray 
illumination, an arrangement which allows us to flt Lx as a free parameter in our model. 

In future observations of the time- variability of the wind lines, the light travel-time 

effects may be used to advantage, as the wind may act as a "low-pass fllter" to the X-ray 
observations, with the fllter cutoff indicating the size of the ionized region (Kallman, McCray, 
& Voit 1987). 
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Table 1: Log of Observations 



Obs. ID: 








Exposure 


Orbital 


HST Archive ID 


Uo. o idi I J- iiiic 


A/TTTt 

iVl J u 






Dataset 1: 


O8HD1010 


2002-06-24 18:10:02 


52449.757 


2676.00 


0.55 


Dataset 2: 


O8HD1020 


2002-06-24 19:36:41 


52449.817 


2172.00 


0.56 


Dataset 3: 


O8HD2010 


2002-06-27 16:39:48 


52452.694 


2172.00 


0.07 


Dataset 4: 


O8HD2020 


2002-06-27 18:06:34 


52452.754 


2676.00 


0.08 


Dataset 5: 


O8NX1010 


2003-07-05 04:27:23 


52825.185 


2147.00 


0.59 


Dataset 6: 


O8NX1020 


2003-07-05 05:51:58 


52825.244 


2675.00 


0.60 


Dataset 7: 


O8NX2010 


2003-07-07 04:27:47 


52827.185 


2147.00 


0.95 


Dataset 8: 


O8NX2020 


2003-07-07 05:53:29 


52827.245 


2675.00 


0.96 



* Phase computed at center of observation. 



Table 2: Line identifications for the spectrum of Figure 4. 





Laboratory 






Line ID 


Wavelength (A) 


Transition 


Comment 


Lya 


1215.670 


1-2 1/2-* 


saturated 


N V 


1238.821 


2q 2po , 
i->- -T 1/2-3/2 


saturated 


N V 


1242.804 


2q 2po , 
ij- -T 1/2-1/2 


saturated 


Si II 


1260.422 


2P°-2D 1/2-3/2 




I 


1302.168 


=^P-3S° 2-1 


weak,blended with Si I I 


I 


1304.858 


3P-3S° i_i 


weak,blendcd with Si II 


Si II 


1304.370 


2po 2q 
r - D 1/2-1/2 


weak,blended with I 


Si II 


1309.276 


2po 2q 
r - D 3/2-1/2 


weak,blended with I 


C II 


1334.532 


2P°-2D 1/2-3/2 


weak 


C II 


1335.663 


2tdo 2rv 
r - U 3/2-3/2 


weak 


Si IV 


1393.755 


2q 2po , 
i->- -T 1/2-3/2 




Si IV 


1402.770 


2c 2po , 
i->- -T 1/2-1/2 




Si II 


1526.707 


2po 2c , 
r - 1/2-1/2 




Si II 


1533.431 


2po 2c , 
-T - S 3/2-1/2 




C IV 


1548.203 


2q 2po 
^- -T 1/2-3/2 




C IV 


1550.777 


2c 2po 
iJ- -T 1/2-1/2 




He II 


1640.420 


2-3 *-* 


weak 


Al II 


1670.787 


^S-ip° 0-1 
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Table 3: Model Pcirameters 

Symbol Adopted value Meaning 



i 40° Orbital inclination 

Ro 1.5 X 10^^ cm Radius of O star 

i?orbit 2.0 Rq Semimajor axis 

Vsini 94 km s~^ (Projected 

rotation velocity) 

Vgys Systemic velocity 

oc 3.6x10"^ Carbon abundance 

On 1.1x10"^ Nitrogen abundance 

osi 3.5x10"^ Silicon abundance 



Comments 

< 20 - 67° (Gies et al. 2003, Tarasov et al. 2003; 
Wen et al. 1999; Herrero et al. 1995; 
Friend & Cassinclli 1986; Davis & Hartmann 1983) 
1.23 X 10^2 cm (Davis & Hartmann 1983) 
1.2 - 1.6 X 10^2 cm (Gies & Bolton 1986) 
2.4 Ro (Davis & Hartmann 1983) 
94.3 ± 5 (Ninkov et al. 1987), 
96 (Gies & Bolton 1986) 
Qtl km s-^ (Ninkov et al. 1987) 
solar from Verner et al. 1994 
solar from Verner et al. 1994 
solar from Verner et al. 1994 



Table 4: Best-fit values for free peirameters for Cygnus X-1 P Cygni line fits 



Symbol Fitted value (Si IV, CIV, NV) 

^ 0.748±0.003 

(0.436±0.007, 0.88±0.03, 0.63±0.05) 

Twind (12.7±0.2, 25.2±0.1, 7.8±1.0) 

(1420±10, 1330±30, 1430) km s"^ 

^;turb (160±10, 292±1, 200) km s"^ 

ai (-0.31±0.03, 0.32±0.01, 0.66±0.10) 

as (0.89±0.06, 0.140±0.004, 0.80±0.01) 

L^^^^/M_Q (0.033±0.001, 0.38±0.01, 0.12±0.01) 

M_6 (4.8±0.3, 4.8±0.1, 6.7±0.2) 

Xl (3.0, 5.8, 1.3) 

V (933, 596, 585) 



Definition* Comments 

Wind acceleration Fit to Si IV, 

then fixed 

Optical depth of stellar 
absorption line 
Wind total depth 
Wind terminal velocity 2300 km s^^ (Davis 

& Hartmann 1983), 
Frozen for N V 
Wind microturbulence Frozen for N V 
Wind opacity exponent 
Wind opacity exponent 
Ratio, X-ray luminosity 
(10^® erg s~^) to wind 
mass loss (10~^ yr~^) 
Wind mass loss rate 
(10-6 yr-i) 
Goodness of fit 
Degrees of freedom 



*From Lamers, Cerruti-Sola, & Perinotto (1987). 
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2600 2800 3000 3200 3400 

JD-2450000.00 



Fig. 1. — One day averages of the flux observed from Cyg X-1 with the AU-Sky Monitor 
on board the Rossi X-Ray Timing Explorer. Arrows indicate the times of the HST/STIS 
observations. 
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1300 1400 1500 1600 1700 



Wavelength (A) 

Fig. 2. — Unsmoothed, raw, HST/STIS data of Cygnus X-1 stacked with increasing orbital 
phase, regardless of epoch. 
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2002 Observations 

= 0.55 0o,b=OO7 



2003 Observations 



8x10-1" 



6x10-1* 



4x10-1'! 



1545 1550 1555 1560 



1540 1545 



1550 




1400 1405 



1555 1560 



Wavelength (A) 



Fig. 3. — STIS observations of Cyg X-1 taken roughly a year apart showing enhanced 
absorption at orbital phases when the X-ray source is behind the companion star (phase 
0.0; solid lines) and illustrating the stability of the absorption change as a function of or- 
bital phase. Narrow features are interstellar. Black dashed lines indicate the laboratory 
wavelengths of transitions identified in Table 2. 
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1300 1400 1500 1600 1700 

Wavelength (A) 



Fig. 4. — A Cyg X-1 STIS spectrum at phase 0.96 with major spectral features indicated. 
The data have been smoothed with a 7-point boxcar function to better show the features of 
interest. Sharp features are interstellar absorption lines. 
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Fig. 5. — Spectra of the N V line at two different orbital phases and two diferent epochs. 
The dashed lines indicate the laboratory wavelengths for transitions listed in Table 2. 
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Fig. 6. — As in figure [5] for Si II 
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Fig. 7. — As in figure [5] for C II 
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Fig. 8. — As in figure [5] for Si IV 
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Fig. 9. — As in figure [5] for C IV 
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Fig. 10. — As in figure [5] for Hell 
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Fig. 11.— Closeup of Cyg X-1 N V line profiles. The blue line depicts a model computed 
using the SEI method (Lamers, Cerruti-Sola, & Perinotto 1987), using a wind velocity law 
given by (1 - 1/r)^. 
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Fig. 12.— As in Fig. [ID for Si IV. 
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Fig. 13.— 



As in Fig. M for C IV. 
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Fig. 14. — The response of the wind-formed P Cygni hnes to X-ray luminosities of l/SL^., 
Lx, and S/SL^. (dotted, bold, and regular thickness, respectively). All other parameters are 
as assumed for the best fit to the Si IV hnes. This figure shows the model at 0=0.55 and 
ignores the effects of hght travel-time, which would tend to decrease the change in the hne 
with changing X-ray luminosity. 



-28- 



2.5 



5/3 Lx 




1390 



1395 
Wavelength (A) 



1400 



1405 



Fig. 15. — The response of the wind- formed P Cygni hnes to X-ray luminosities of l/SL^., 
Lx, and S/SL^. (dotted, bold, and regular thickness, respectively). All other parameters are 
as assumed for the best fit to the Si IV hnes. This figure shows the model at 0=0.96 and 
ignores the effects of hght travel-time, which would tend to decrease the change in the hne 
with changing X-ray luminosity. 
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Fig. 16. — Contours of Logio of Si IV abundance outside of the shadow region. We assume 
the X-ray luminosity is (1/3)1.6 x 10^^ erg s^^, and that all other parameters have their best- 
fit values. We include only ionization from the black hole and not the ambient ionization of 
the wind. The circle represents the supergiant and the X the position of the compact source. 
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Fig. 17. — Contours of Logio of Si IV abundance outside of the shadow region. We assume 
the X-ray luminosity is 1.6 x 10^'' erg s~^, and that all other parameters have their best-fit 
values. We include only ionization from the black hole and not the ambient ionization of the 
wind. 



X 



Fig. 18. — Contours of Logio of Si IV abundance outside of the shadow region. We assume 
the X-ray luminosity is (5/3)1.6 x 10"^^ erg s^^, and that all other parameters have their best- 
fit values. We include only ionization from the black hole and not the ambient ionization of 
the wind. 
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Fig. 19. — Contours of radial optical depth r, given the X-ray luminosity is (1/3)1.6 x 
10^''' erg 
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Fig. 21. — Contours of radial optical depth r, given the X-ray luminosity is (1/3)1.6 x 
10^''' erg 



